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ABSTRACT 

We adopt a new chemical evolution model for the Large Magellanic Cloud (LMC) and thereby in- 
vestigate its past star formation and chemical enrichment histories. The delay time distribution of 
type la supernovae recently revealed by type la supernova surveys is incorporated self-consistently into 
the new model. The principle results are summarized as follows. The present gas mass fraction and 
stellar metallicity as well as the higher [Ba/Fe] in metal-poor stars at [Fe/H]< —1.5 can be more self- 
consistently explained by models with steeper initial mass functions. The observed higher [Mg/Fe] 
(> 0.3) at [Fe/R] 0.6 and higher [Ba/Fe] (> 0.5) at [Fe/H] 0.3 can be due to significantly en- 
hanced star formation about 2 Gyr ago. The observed overall [Ca/Fe]-[Fe/H] relation and remarkably 
low [Ca/Fe] (< —0.2) at [Fe/H]> —0.6 are consistent with models with short-delay supernova la and 
with the more efficient loss of Ca possibly caused by an explosion mechanism of type II supernovae. 
Although the metallicity distribution functions do not show double peaks in the models with a starburst 
about 2 Gyr ago, they show characteristic double peaks in the models with double starbursts at ~ 200 
Myr and ^ 2 Gyr ago. The observed apparent dip of [Fe/H] around ~ 1.5 Gyr ago in the age-metallicity 
relation can be reproduced by models in which a large amount (~ IO^Mq) of metal-poor ([Fe/H]< —1) 
gas can be accreted onto the LMC. 

Subject headings: Magellanic Clouds ~ stars: abundances - galaxies: evolution 



1. INTRODUCTION 

The time evolution of the chemical abundances in the 
interstellar medium (ISM), field stars, and globular clus- 
ters (GCs) of the Large Magellanic Cloud (LMC) contains 
valuable information on its long-term star formation his- 
tory and thus has been investigated in detail by many au- 
thors (e.g., Da Costa 1991; Olszewski et al. 1991; Russel 
& Dopita 1992; Dopita et al. 1997; Geisler et al. 1997; 
Pagel & Tautvaisiene 1998, PT98; Cole et al. 2005, C05; 
Cioni et al. 2006). The elemental abundance ratios of a 
{a means alpha-elements), Fe-peak, and neutron-capture 
elements in field stars and GCs of the LMC have been ex- 
tensively investigated by spectroscopic observations. The 
observed differences in the abundance ratios between the 
LMC and the Galaxy have been discussed in detail (e.g.. 
Hill et al. 1995, 2000; Johnson et al. 2006). The radial 
and azimuthal variations of stellar metallicities in the the 
LMC disk has been investigated both observationally and 
theoretically in terms of its star formation and dynamical 
evolution histories (Geisler et al. 2003; Bekki & Chiba 
2005; Cioni et al. 2006). 

One of the importance results in these previous stud- 
ies is that the observed age-metallicity relation (AMR) is 
more consistent with a model with a secondary "starburst" 
about a few Gyr ago in the LMC (e.g., Tsujimoto et al. 
1995, T95; PT98). Although the possible presence of a 
past starburst (or significantly enhanced star formation) 
has long been discussed in other observational studies (e.g.. 
Butcher 1977; Stryker 1983; Bica et al. 1986; Bcrtelh et 
al. 1992; Olszewski et al. 1996; Gallagher et al. 1996; 



Vallenari et al. 1996; Ardeberg et al. 1997; Elson et al. 
1997; Geha et al. 1998 Holzman et al. 1998, 1999; Olsen 
1999; Smecker-Hane et al. 2002; Glatt et al. 2010; Indu 
& Subramaniam 2011), the epoch and the strength of the 
burst were not well constrained. Furthermore, owing to 
the lack of modern chemical evolution models with the 
latest chemical yields of asymptotic giant branch (AGB) 
stars and supernovae, it remained unclear how the long- 
term star formation history with a possible starburst could 
be imprinted on the detailed elemental abundance ratios 
of field stars and GCs in the LMC. 

Recent photometric and spectroscopic observations of 
field stars and GCs have provided new clues for these un- 
resolved problems in the LMC. For example, photometric 
studies of the stellar populations in the LMC have revealed 
the AMRs of different local regions in the LMC disk (e.g., 
Harris & Zaritsky 2009, HZ09; Rubele et al. 2012, R12). 
The AMR derived by HZ09 show enhanced star formation 
rates of the LMC around 2 Gyr, 500 Myr, 100 Myr, and 
12 Myr with the two peaks (500 Myr and 2 Gyr) being 
nearly coincident with the star formation peaks observed 
in the SMC. Piatti (2011) has also shown that there was a 
burst of cluster formation around 2 Gyr in the LMC and 
suggested that the LMC experienced strong tidal inter- 
action with the SMC and possibly with the Galaxy. Us- 
ing deep near-infrared data from the VISTA near-infrared 
YJKg survey of the Magellanic system, R12 derived the 
star formation histories of different local regions in the 
LMC disk. They revealed the presence of peaks in SFRs 
around 3 Gyr and 5 Gyr ago for most of the subregions. 
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Glatt et al. (2010) have recently found an enhanced clus- 
ter formation at 125 Myr and 800 Myr ago in their 324 
populous star clusters for the LMC. 

Furthermore, recent spectroscopic observations of field 
stars and GCs in the LMC have revealed intriguing chem- 
ical properties of the LMC (e.g., Mucciarelli et al. 2008, 
Pompeia et al. 2008, P08; Colucci et al. 2012, C12, 
Haschke et al. 2012). Mucciarelli et al. (2008) have shown 
that all four of the intermediate-age GCs which they in- 
vestigated in the LMC have negligible star-to-star scatter 
in their chemical abundances of light, a, iron-peak and 
neutron-capture elements. This implies that secondary 
star formation from gaseous cjccta of stars within the GCs 
did not occur so that the chemical abundances might be 
similar between intermediate- age field stars and clusters. 
P08 have found that [Ca/Fe] and [Si/Fe] in the LMC stars 
are lower than those of the Galactic stars at the same 
[Fc/H] whereas the [0/Fe] and [Mg/Fc] arc only slightly 
deficient in comparison with their Galactic counterparts. 
C12 have revealed the age-dependence of [a/Fo] and [Fe/H] 
among globular /star clusters with ages between 0.05 and 
12 Gyr and the significant enhancement of the neutron- 
capture elements Ba, La, Nd, Sm, and Eu in the youngest 
star clusters. The origin of these recent observational re- 
sults has not yet been treated by chemical evolution mod- 
els. 

In spite of this progress in observational studies of the 
stellar populations of the LMC, theoretical models to ex- 
plain the observations have not yet been fully developed. 
T95 tried to explain the observed [0/Fe]-[Fe/H] relation, 
the metallicity distribution function (MDF) , and the AMR 
in a self-consistent manner by adopting a model in which 
the LMC had a starburst about 3 Gyr ago and a steeper 
initial mass function (IMF). They explained these obser- 
vations, though the observational data points are much 
smaller than those that we can now access. Pilyugin (1996) 
demonstrated that the observed [Fe/0] versus [Fe/H] re- 
lation can be reproduced only by the galactic wind model 
in which stellar ejecta can be preferentially expelled from 
the LMC. PT98 explained the observed AMR and chemi- 
cal abundance patterns of the LMC by considering "non- 
selective stellar wind" models in which some fraction of 
the stellar ejecta both from AGB stars and SNe can be ex- 
pelled completely from the LMC and consequently can not 
participate in the chemical enrichment processes. These 
previous chemical evolution models did not incorporate 
the delay time distribution (DTD) of type la supernova 
(SNe la) recently revealed by extensive SN la surveys (e.g., 
Mannucci et al. 2006; Sullivan et al. 2006; Totani et al. 
2008; Maoz et al. 2010) and thus might not be regarded as 
realistic and reasonable. The construction of more sophis- 
ticated and realistic models is indispensable for discussing 
the above mentioned latest observational results. 

The purpose of this paper is to adopt a new chemical 
evolution model and thereby discuss the latest observa- 
tional results on the chemical evolution and star forma- 
tion history of the LMC. The observed DTD of SNe la is 
adopted in the new model so that the progenitor stars of 
SNe la can explode as early as 10**yr after their formation. 
These "prompt SNe la" can cause fundamental differences 
in the chemical evolution between the present model and 
previous ones (e.g., T95) in which the time delay between 



star formation and explosion of SNe la is assumed to be 
typically ~ 10^ yr ("classical SNe la"). The new model 
also incorporates the metallicity dependent chemical yields 
of AGB stars (e.g., Busso et al. 2001; Tsujimoto & Bekki 
2012; TB12) so that the chemical evolution of s-process 
elements can be investigated and compared with the ob- 
servations. The new model is the so-called "one-zone" 
model, in which the time evolution of averaged chemical 
properties can be investigated. 

Wc; explore a wide range of models with and without 
starburst and different star formation histories in the LMC 
so that we can discuss recent new observational results on 
different chemical properties of the LMC disk stars. For 
example, we will discuss (i) whether and how the previous 
starburst (s) triggered by the LMC-SMC-Galaxy interac- 
tion can be imprinted on the chemical abundance proper- 
ties of the LMC, (ii) how the IMF can control the chem- 
ical evolution, and (iii) how the accretion of metal-poor 
gas from the SMC or high velocity clouds onto the LMC 
(Bekki & Chiba 2007, BC07; Diaz & Bekki 2012, DB12; 
Bekki & Tsujimoto 2010) influences the chemical evolu- 
tion. It is timely to discuss the chemical evolution of the 
LMC in the context of the Galaxy LMC SMC tidal in- 
teraction, given that recent observational and theoretical 
studies have demonstrated that the tidal interaction is im- 
portant not only for the star formation history of the LMC 
but also for the stellar and gaseous distributions of the 
Magellanic system (e.g., Olsen et al. 2011; DB12). 

The layout of this paper is as follows. In §2, we describe 
our new one-zone chemical evolution models in detail. In 
§3, we present the results of the time evolution of the chem- 
ical abundances for models with different parameters. In 
§4, we discuss the results in the context of the IMF, the 
long-term star formation history, and the gas accretion 
events in the LMC. The conclusions of the present study 
are given in §5. 

2. THE MODEL 

2.1. Basic equations 

We adopt one-zone chemical evolution models that are 
essentially the same as those adopted in our previous stud- 
ies on the chemical evolution of the Galaxy (Tsujimoto 
et al. 2009; TB12). Accordingly, we briefly describe the 
adopted models in the present study. The present model is 
improved in comparison with previous ones (Tsujimoto et 
al. 2009) in terms of including the s-process elements and 
prompt SNIa self-consistently in the models. The LMC 
disk is assumed to form through a continuous gas infall 
from outside the disk region (e.g., halo) for the last 13 
Gyr, as in previous chemical evolution models (e.g., Tsu- 
jimoto et al. 2009). Some fraction of gas with metals can 
be expelled from the LMC as stellar winds through ener- 
getic feedback effects of SNe in some models of the present 
study (e.g., "wind models", as described later). The star 
formation can be suddenly and significantly enhanced (re- 
ferred to as "bursts of star formation" or "starburst"). As 
demonstrated by theoretical studies, the Galaxy-LMC- 
SMC tidal interaction is important both for the formation 
of the Magellanic Stream and for reproducing the observed 
peaks of star formation in the LMC (e.g., Bekki & Chiba 
2005; DB12). Therefore, the adopted assumption of star- 
burst is quite reasonable. Although the infall of metal- 
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poor gas from the SMC and other dwarf galaxies can be 
important for the LMC chemical evolution, we discuss this 
later in §4. 

We investigate the time evolution of the gas mass frac- 
tion (/g(t)), the star formation rate {'ip{t)), and the abun- 
dance of the ith heavy element {Zi{t)) for a given accretion 
rate {A{t)), IMF, and ejection rate of ISM {w{t)). The ba- 
sic equations for the adopted one-zone chemical evolution 
models are described as follows: 

^^-aMt)+A{t)-w{t) (1) 

= + ZA,i{t)A{t) + yu,im 

+yu,i / tp{t - tu)g{tu)dtia 
Jo 

+ / 2/agb,i(TOagb)V'(i-iagb)^(tagb)(^iagb - Wi(i) , (2) 

Jo 

where ai is the mass fraction locked up in dead stellar 

remnants and long-lived stars, j/ia,i, 2/ii.i, and j/agb.i are the 
chemical yields for the ith clement from type II supernovae 
(SN II), from SN la, and from AGB stars, respectively, 
ZA,i is the abundance of heavy elements contained in the 
infalling gas, and Wj is the wind rate for each element. The 
quantities tia and tagb represent the time delay between 
star formation and SN la explosion and that between star 
formation and the onset of AGB phase, respectively. The 
terms fl'(iia) and h^gh are the distribution functions of SNe 
la and AGB stars, respectively, and the details of which 
are described later in this section. The term ftagb con- 
trols how much AGB ejecta can be returned into the ISM 
per unit mass for a given time in equation (2). The to- 
tal gas masses ejected from AGB stars depends on the 
original masses of the AGB stars (e.g., Weidemann 2000). 
Therefore, this term /lagb depends on the adopted IMF and 
the age-mass relation of the stars (later described). Thus 
equation (1) describes the time evolution of the gas due to 
star formation, gas accretion, and stellar wind. Equation 
2 describes the time evolution of the chemical abundances 
due to chemical enrichment by supernovae and AGB stars. 

The star formation rate tp{t) is assumed to be propor- 
tional to the gas fraction with a constant star formation 
coefficient and thus is described as follows: 

m = cMt) (3) 

We assume that Cgf is different between (i) the "quies- 
cent phase" , when the LMC shows an almost steady star 
formation, and (ii) the "burst phase" , when the LMC ex- 
perienced a burst of star formation. In the present study, 
we investigate models with no burst (referred to as "stan- 
dard models" ) and those in which a starburst can occur 
only once ( "burst models" ) and twice ( "double-burst mod- 
els"). The star formation rate is assumed to be constantly 
higher for fgbi.s < t < fgbi.e in the first starburst and 
^sb2.s < t < tsb2,c in the second one. The star formation 
coefficient (Cgf) is thus described as follows: 

{Cq for quiescent phase 
Csbi for ffist starburst (4) 
Csb2 for second starburst 

In the present study, a non-dimensional value of Cgf is 
given for each model. 



For the accretion rate, we adopt the formula in which 
A{t) = Caexp(— i/ia) and tg, is a free parameter control- 
ling the time scale of the gas accretion. The normalization 
factor Ca is determined such that the total gas mass ac- 
creted onto the LMC can be 1 for a given ta. Although we 
investigated models with different ts,, we finally adopt the 
models with tg, = 0.3 Gyr as reasonable ones for the LMC. 
This is mainly because the models with ta = 0.3 Gyr can 
better explain the observations. The present results do not 
depend strongly on the parameter ia for a reasonable pa- 
rameter range. The initial [Fe/H] of the infalling gas is set 
to be —2 and we assume a SN-II like enhanced [a/Fe] ratio 
(e.g., [Mg/Fe]w 0.4) and [Ba/Fe] for the gas. This adop- 
tion of initial [Fe/H] is reasonable for the present study 
which discusses stars with [Fe/H] as low as —2. Also, such 
adoption has been done in other chemical evolution models 
(e.g., Tsujimoto et al. 2009). We discuss our parameter 
study for the above variables in §2.5. 

2.2. Selective and non-selective stellar winds 

We investigate models in which metals from AGB stars 
and SNe can be removed from the LMC through energetic 
stellar winds so that they can not be used for chemical 
enrichment of the LMC ISM. These models are referred to 
as "wind models" for convenience. The wind models are 
further divided into two categories: "selective wind mod- 
els" and "non-selective wind models". Only SNe ejecta 
(not already existing ISM) can be expelled from the LMC 
in the selective wind model, whereas both ejecta from SNe 
and AGB stars and already existing ISM can be expelled 
from the LMC in the non-selective wind model. For com- 
parison, we also investigate models in which there is no 
stellar wind and thus gaseous ejecta from all stars can be 
fully mixed with the ISM ( "non-wind models" ) . 

In the selective wind models, some fraction, 1 — /ej, of 
gaseous ejecta from SNe can be mixed with ISM for chem- 
ical enrichment processes whereas all AGB ejecta can be 
mixed with ISM. Therefore, the wind rate w(t) is esti- 
mated only from the total mass of gaseous ejecta from 
SNe Mej^sn at each time step. The total mass, M^, of ISM 
that is removed from the LMC at each time step in the 
selective wind models is 

Mw = /ejMej,,n. (5) 

In the non-selective wind models, we adopt the same 
model as that used in PT98 in which Mej is determined 
solely by the star formation rate: 

Mw = CejV(t), (6) 

where Cej is a parameter (in simulation units) that can 
control how efficiently the ISM with AGB and SN ejecta 

can be removed from the LMC. We investigate non-wind 
(i.e., standard, burst, non-burst), selective wind, and non- 
selective wind models. We discuss our parameter study 
for the above variables in §2.5. 

2.3. Chemical yields and delay time distribution of SN la 

We adopt the nucleosynthesis yields of SNe II and la 
from T95 to deduce yu,i and yia,i for a given IMF. Stars 
with masses larger than 8Mq explode as SNe II soon af- 
ter their formation and eject their metals into the ISM. In 
contrast, there is a time delay (iia) between the star for- 
mation and the metal ejection for SNe la. We here adopt 
the following DTD (^(iia)) for 0.1 Gyr < tu < 10 Gyr, 
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which is consistent with recent observational studies on 
the SN la rate in extra-galaxies (e.g., Totani et al. 2008; 
Maoz et al. 2010, 2011): 

5la(tla) = (7) 

where Cg is a normalization constant that is determined 
by the number of SN la per unit mass (which is controlled 
by the IMF and the binary fraction for intermediate-mass 
stars for the adopted power-law slope of —1). Maoz & 
Badens (2010) detected a population of prompt SNIa in 
the LMC and showed that the number of prompt SNIa per 
stellar mass formed is 2.7— 11.0 x IO^^Mq^. Although we 
mainly investigate the "prompt SN la models" with the 
above DTD, we also investigate the "classical SN la mod- 
els" with 1 Gyr < tu < 3 Gyr (e.g., Yoshii et al. 1996). 
This SNIa lifetime in Yoshii et al. (1996) was deduced by 
using their Galactic chemical evolution models that can 
be consistent with the observed [0/Fe]-[Fe/H] relations of 
the Galactic stars in the solar neighborhood. The fraction 
of the stars that eventually produce SNc la for 3-8Af0 has 
not been observationally determined and thus is regarded 
as a free parameter, /b. We mainly investigate models 
with /b = 0.03, because such models can better explain 
the observed chemical properties of the LMC. We briefly 
compare the results of models with /b = 0.03 and 0.06. For 
the site of r-process, we adopt the mass range of S-IOM© 
for SNe II (Mathews et al. 1992; Ishimaru et al. 1999) as 
identified with the collapsing 0-Ne-Mg core (Wheeler et 
al. 1998). The yield of Ba from r-process is 1.45 x IO^^Mq 
for the adopted range of SNe II. 

Low-mass AGB stars (< 3Mq) release the s-process 
elements during the thermally pulsing AGB phase (e.g., 
Gallino ot al. 1998). As a consequence of large uncer- 
tainties in convectivc mixing and ^"^C-pocket efficiencies, 
the s-process nucleosynthesis allows for a wide range of 
possible production levels. On the other hand, the ob- 
served abundances for the AGB stars can be directly com- 
pared with the theoretical nucleosynthesis results. Here 
we investigate the element Ba by adopting the best empir- 
ical metallicity dependent Ba-yield derived by TB12. The 
adopted metallicity dependent Ba-yield j/Ba is 

„ _ / 2/Ba,o(magb) X lO^MFe/U]+1.5 _2 < [Pe/H] < -1 

^""^ 1 yBa,o(magb) -1 < [Fe/H] 

(8) 

The value of T/Ba,o("i'agb) depends on the mass of the 
stars magb that can finally become AGB stars. The 
adopted yBa,o('Tiagb) are 2.6 x 10~'^Mq, 4.0 x IQ-'^M©, 
and 5.2 x IQ-'^Mq for rriagb = l-SM©, 2.OM0, and 3.0Mq, 
respectively. 

An AGB star with initial mass mj and final mass m-p 
can eject its envelope with a total mass of rriej and the 
gas can be mixed with the surrounding ISM to chemically 
pollute the ISM. The initial-to-final relationship for AGB 
stars, from which we can deduce mej , has been extensively 
discussed in Weidemann (2000). We derive an analytic 
form for mej (= mi — mp) from the observational data 
by Weidemann (2000) by using the least squares fitting 
method, and find 

mej = O.9I6M1 - 0.444. (9) 

The coefficient of determination (i?-value) is 0.995 for the 
above fitting. In order to calculate the main-sequence 



turn-off mass mxo ; we use the following formula (Renzini 
& Buzzoni 1986): 

logmTo(ts) = 0.0558(logis)^ - 1.338 logts -1-7.764, (10) 

where mxo is in solar units and time is in years. By 
using the adopted IMF and equations (9) and (10), we nu- 
merically estimate h^gh (i-e., how much AGB ejecta can be 
returned into ISM per unit mass) in equation (2) at each 
time step. 

2.4. IMF 

The adopted IMF is defined as ^'(toi) = Ms,o"^I"", 
where mi is the initial mass of each individual star and the 
slope a = 2.35 corresponds to the Salpeter IMF (Salpeter 
1955). The normalization factor Ms ft is a function of a, mi 
(lower mass cut-off), and (upper mass cut-off). These 
mi and are set to be O.IMq and 50Mq, respectively 
(so that the normalization factor M^.o is dependent simply 
on a). We investigate models with different a to find the 
model(s) that can best explain the observed abundance 
patterns of stars in the LMC. We do not discuss models 
with different m^, because the effects of changing m^ on 
the LMC chemical evolution are similar to those of chang- 
ing a. 

2.5. Parameter studies 

Wc mainly investigate the standard (labeled as "S"), 
burst ("B"), and double-burst ("DB") models in which 
stellar wind effects are not included at all (i.e., "non-wind 
models"). By using these non-wind models, we demon- 
strate how the IMF slope (a) and the epochs of starburst 
(fsbi.s and isb2.s) can infiucncc the chemical evolution of 
the LMC. These tsbi.s and tsb2,s are given in units of Gyr. 
We then investigate the wind models ("W") so that we 
can discuss whether or not removal of AGB and SN ejecta 
from the LMC is important in the chemical evolution of 
the LMC. In the present study, the time t is the time that 
has elapsed since the model calculation started. Therefore 
t = Gyr (13 Gyr) represent the time when the calcula- 
tion starts (ends). Previous observational and theoretical 
studies have suggested that there could be at least two 
epochs of enhanced star formation (starburst) about ^ 2 
Gyr ago and 0.2 Gyr ago (e.g., HZ09 and DB12). We 
therefore mainly investigate the burst and double-burst 
models with fgbi.s = 11 Gyr and fsb2.s = 12.8 Gyr. Table 
1 summarizes the representative 26 models investigated in 
the present study. 

Figure 1 illustrates the time evolution of star formation 
rates (SFRs) and [Fe/H] of stars in the five representative 
models (SI, Bl, B6, DBl, and DB6) with different pa- 
rameters controlling the LMC star formation history. No 
burst of star formation is assumed in SI whereas only one 
burst is assumed in Bl and B6. Two bursts of star for- 
mation at different epochs are assumed in DBl and DB6. 
These models are chosen so that they can, in combination, 
show a wide range of star formation histories in Figure 1: 
they are just examples of possible star formation histories 
of the LMC. We investigate how the final chemical abun- 
dances depend on the LMC star formation history by using 
the results of models with different star formation histo- 
ries. Figure 2 shows the time evolution of [Mg/Fe] and the 
[Mg/Fe]-[Fe/H] relation in the models with prompt and 
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classical SNe la. It is clear from this figure that the models 
with prompt SN la show no plateau in the [Mg/Fe]-[Fe/H] 
relation owing to the earlier chemical enrichment by SNe 
la that can eject an Fe-rich gas. Figure 2 also shows that 
the time evolution of [Mg/Fc] depends on /b in such a way 
that [Mg/Fe] can steeply decrease with time in the model 
with larger /b- We discuss these in more detail for different 
models in the following sections. 

2.6. Observations to be compared with predictions 

We mainly investigate (i) the present gas mass fraction 
/g,o, (ii) the stellar metallicity of the youngest stellar pop- 
ulation, [Fe/H]o, (iii) AMR, (iv) [Mg/Fc]- [Fc/H] relation 
(as an example of [Q;/Fe]-[Fe/H] relations), (v) [Ba/Fe]- 
[Fe/H] relation (as an example of the time evolution of 
s— process elements), and (vi) MDF. We compare these re- 
sults with the corresponding observations for the LMC. In 
order to demonstrate clearly how different chemical prop- 
erties of the LMC are compared with those of the Galaxy, 
we also show the observational results for the Galaxy. We 
estimate the total stellar mass (Mg) of the LMC by using 
the observed F-band luminosity (« 3 x 10^ Lq) and the 
reasonable mass-to-light ratio of 0.9 ± 0.2 for the observed 
B — V color (Bell & de Jong 2001). By using the observed 
total mass (Mg) of the LMC ISM (Bernard et al. 2008) 
and Ms, we can estimate /g,o- The observational error 
bar of /g,o is due largely to the uncertainty of the stellar 
mass-to-light ratio for the LMC. We adopt [Fe/H]o ~ —0.3 
as derived by Luck et al. (1998) for Cepheid variables 
with ages of 10-60 Myr, because these youngest popula- 
tions can have the present-day stellar metallicity of the 
LMC. Each of the observed AMRs shows a wide spread 
in each age bin (represented by an error bar) owing to 
the presence of stars with different [Fe/H] at each age bin. 
Although this could make it difficult for us to derive the 
best model, wc try to give at least some constraints on 
some of the model parameters. Wc mainly investigate time 
evolution of [Mg/Fe] and [Ba/Fe], because these abun- 
dances are more reliably derived from observations. The 
[Mg/Fe]-[Fe/H] and [Ba/Fe]-[Fe/H] relations are examples 
of [a/Fe] [Fe/H] and [s-process/Fe] relations, respectively, 
in this study. 

3. RESULTS 

3.1. Standard models 
3.1.1. AMR 

First we describe the five standard models with different 
IMFs in order to demonstrate the importance of the IMF 
slopes in the LMC chemical evolution. In these models, Cq 
is chosen such that the final stellar metallicity ([Fe/II]o) 
can be —0.3 (i.e., the observed value). Figure 3 shows the 
AMRs for the five models as well as the observed [Fe/H] 
of the LMC field stars at different age bins and and in- 
dividual GCs. The AMR in each model is simply a plot 
of [Fe/H] at each time step (i.e., not average [Fc/H] over 
a given age bin like observation). Although the observa- 
tional errors are not small for old GCs in the LMC, their 
data points are included to compare stars with the lowest 
[Fe/H] in the models with the corresponding observations. 

The five models with no burst of star formation can 
reproduce reasonably well the overall trend of the ob- 
served AMR, which implies that the AMR alone can not 



be used for discriminating burst and non-burst models for 
the LMC. These models, however, appear to be less con- 
sistent with the observed [Fe/H] around ages of 2-3 Gyr 
(i.e., t =10-11 Gyr) owing to the presence of stars with sig- 
nificantly lower metallicitics (—1 < [Fc/H] < —0.5). The 
models with shallower IMFs show higher [Fe/H] at a given 
age, which reflects the fact that chemical evolution can 
proceed more rapidly owing to a larger amount of metals 
produced by a larger number of SNe. 

It shoidd bo stressed that these standard models can 
not reproduce so well the AMR around ages of 3-9 Gyr 
(i.e., t =4-10 Gyr) derived by HZ09 in which the AMR 
shows systematically lower [Fc/H] for a given age in com- 
parison with other observations by C05 and Carrera et al. 
(2008, C08). The reason for this apparent difference in the 
observed AMRs between different observations could be 
related to different methods to determine ages and metal- 
licitics of stars in the observations. If the results by HZ09 
are closer to the true AMR of the LMC, then the stan- 
dard models with no burst can be regarded as less realistic 
models for the LMC evolution. Although the large [Fe/H] 
dispersion around an age of 2 Gyr (i.e., t = 11 Gyr) could 
be simply due to star formation from gas with different 
metallicities (i.e., inner higher and outer lower metallici- 
ties of the LMC ISM), it could also be caused by a sudden 
and rapid infall of metal-poor gas from outside the LMC 
disk. If the observed dispersion is due to an external gas 
infall, then it would have profound implications for the 
LMC evolution. We will later discuss the implications in 
§4. 

3.1.2. [Mg/Fe] 

Figure 4 shows that [Mg/Fe] slowly and monotonically 
decreases with time regardless of the adopted IMFs. This 
time-dependence of [Mg/Fe] is simply a result of a later 
contribution of SNIa to the chemical evolution of the LMC. 
The more significant decrease of [Mg/Fc] with time can 
be seen in the models with steeper IMFs owing to the 
greater contribution of SNIa to the chemical evolution in 
these models. The [Mg/Fe] -[Fe/H] relations do not show 
a plateau at low [Fe/H]. The lack of a plateau in the 
[Mg/Fe] [Fe/H] relations reflects the fact that SNe la can 
chemically pollute the LMC ISM from a very early stage 
{t ~ 10^ yr) of its evolution. The lack of a plateau ap- 
pears to be seen in the observed [Mg/Fe]-[Fe/H] relation 
for [Fe/H]< -1.5, though it is not so clear. The [Mg/Fe]- 
[Fe/H] relations in the models with a < 2.75 can reproduce 
reasonably well the overall trend of the observed [Mg/Fe] 
relation, though the observation shows a large [Mg/Fe] 
dispersion at a given [Fe/H]. 

The model with a very steep IMF (a = 2.95) shows 
systematically lower [Mg/Fe] and therefore can not repro- 
duce so well the locations of stars in the [Mg/Fe] [Fe/H] 
relation. The models with shallower IMFs show system- 
atically larger [Mg/Fe] in the present standard models. 
Owing to the observed large dispersion of [Mg/Fe], it is 
currently difficult to determine which IMFs {a = 2.35 or 
2.55 or 2.75) can better explain the observed [Mg/Fe]- 
[Fe/H] relation. The observed high [Mg/Fe] (> 0.4) for 
higher [Fe/H] (> —0.7) can not be explained by any of the 
standard models in the present study. 
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3.1.3. [Ba/FeJ 

Figure 5 shows that [Ba/Fc] starts to slowly increase 
with time about 2 Gyr after the commencement of active 
star formation in the LMC disk. This slow [Ba/Fe] in- 
crease is due to the ejection of s— process elements from 
low- mass (< SMq) AGB stars. The observed systemati- 
cally high [Ba/Fe] (> for most stars) at [Fe/H] < -1.5 
seems to be better reproduced by the models with steeper 
IMFs, though the very high [Ba/Fe] (> 0.2) in some stars 
at such low metallicities can not be reproduced at all by 
any of the standard models in the present study. The 
reason for the higher [Ba/Fe] at [Fe/H]< —1 in steeper 
IMFs is that the numerical ratio of SNe II with masses of 
~ 8Mq to those with ~ lOM© (and thus the contribution 
of these SNe to the chemical evolution) can change in the 
steeper IMF and consequently [Ba/Fe] can increase more 
significantly in the early history of the chemical evolution 
of the LMC. The maximum [Ba/Fe] in the model with 
the Salpeter IMF is at most [Ba/Fe]~ 0.4, whereas most 
of the LMC stars at [Fe/H]> -0.6 show [Ba/Fe]> 0.4. 
These results strongly suggest that the LMC could have 
had a steeper IMF (at least steeper than the Salpeter IMF 
with a = 2.35) in its star formation history (if the LMC 
has not experienced starbursts). The observed very high 
[Ba/Fe] (> 0.8) of some stars at [Fe/II]~ —0.6 can not 
be explained at all by any of the standard models in the 
present study. 

3.1.4. fg-[Fe/H] relation 

In the present study, the Cq for each model is cho- 
sen such that the final [Fe/H] (i.e., [Fe/H]o) is consistent 
with the observed one (—0.3). A smaller amoimt of met- 
als can be ejected from SNe in the models with steeper 
IMFs owing to the smaller numbers of SNe. Therefore, 
ISM in the models with steeper IMFs needs to be more 
rapidly consumed by star formation and chemically en- 
riched by SN ejecta to reach [Fe/H]o = —0.3. Figure 6 
shows that (i) the models with steeper IMFs show smaller 
fgfi and (ii) the model with a = 2.55 is the most consis- 
tent with the observed /g.o- These results suggest that the 
observed /g,o and [Fe/II]o combine to support the steeper 
IMF {a. « 2.55) of the LMC. However this suggestion de- 
pends on the adopted assumption that no ISM can be ex- 
pelled from the LMC in the standard models. If a signifi- 
cant fraction of cold ISM can be stripped from the LMC, 
even the models with shallower IMFs (e.g., a = 2.15) may 
be able to explain both /g.o and [Fe/H]o. Likewise, if a 
significant fraction of cold ISM can be recently accreted 
onto the LMC, the models with rather steep IMFs (e.g., 
a > 2.75) may explain both /g^o and [Fe/H]o too. 

3.2. Burst models 
3.2.1. AMR 

The burst of star formation and the subsequent efficient 
production of metals can rapidly and significantly increase 
[Fe/H] in the burst models. Therefore, the chemical evolu- 
tion in the quiescent phase of star formation in the burst 
models needs to proceed more slowly (i.e., smaller Cq) 
in comparison with the standard models so that the final 
[Fe/H] can be as low as —0.3. Figure 7 shows the AMRs 
in the five burst models with difl[erent IMFs and Cq. The 



burst model B3 can not reproduce the observed AMR ow- 
ing to the systematically low [Fe/H] in the quiescent phase 
of star formation. The model Bl with the Salpeter IMF 
can better explain the AMR of C05 and COS whereas the 
models with steeper IMFs (B2 and B5) can explain bet- 
ter the AMR by HZ09. The observed AMR thus cannot 
give strong constraints on the IMF of the LMC. Owing to 
the observed large [Fe/H] dispersion at an age of ~ 2 Gyr 
(t = 11 Gyr), the AMR alone can not allow us to make 
a robust conclusion as to whether the LMC experienced 
a burst of star formation at an age of 2 Gyr (i.e., t = 11 
Gyr). 

3.2.2. [Mg/Fe] 

As shown in Figure 8, the overall trends of the [Mg/Fe] 

evolution in the quiescent phase for the burst models are 
similar to those for the standard ones. The [Mg/Fe]- 
[Fe/H] relations in the burst models show "bumps" (a 
sharp increase followed by a decrease in the time evolu- 
tion of [Mg/Fe]) whose magnitudes depend on Csi and a. 
The observed locations of the LMC stars in the [Mg/Fe]- 
[Fe/H] plane do not show clearly such a bump. However, 
if the two stars with [Mg/Fc]-'o.5 at [Fe/H]~ -1.2 and 
— 1.0 are removed from Figure 8, then the remaining stars 
appear to show a bump around —0.7 <[Fe/H]< —0.6. Fur- 
thermore, the apparent lack of stars with [Mg/Fe] > 0.3 
at [Fe/H]~ —0.4 suggests that [Mg/Fe] has been decreas- 
ing since [Fe/H]= -0.5. Thus the higher [Mg/Fe] at 
[Fe/H]~ —0.6 could be evidence of a star burst around 
t = 11 Gyr (i.e., 2 Gyr ago). However, the observed 
[Mg/Fe] dispersion at [Fe/H]^ —0.6 is so large that we 
can not make a robust conclusion as to whether a star- 
burst occurred in the LMC about t = 11 Gyr (i.e., 2 Gyr 
ago). Also it should be noted that the smaller number of 
stars at [Fe/H]> —0.5 could be responsible for the appar- 
ent lack of stars with higher (> 0.3) [Mg/Fe]. 

3.2.3. [Ba/FeJ 

Figure 9 shows that irrespective of the model param- 
eters, [Ba/Fe] rapidly decreases soon after the starbursts 
owing to the ejection of Fe-rich gas from prompt SNe la 
in the time evolution of [Ba/Fe]. The temporary [Ba/Fe] 
decrease is subsequently followed by its rapid and sharp 
increase due to chemical pollution by AGB ejecta. Ow- 
ing to the rapid [Ba/Fe] increase, the burst models can 
show systematically higher final [Ba/Fe] in comparison 
with the standard models with no burst. Therefore, the 
observed presence of stars with higher [Ba/Fe] (> 0.5) 
at [Fe/H]R:! —0.3 could be more consistent with the burst 
models. The burst models with steeper IMFs (a = 2.75) 
show [Ba/Fe]~ 0.7 at [Fe/H]= -0.3, which suggests that 
the combination of a steeper IMF and a secondary star- 
burst could be closely associated with the origin of stars 
with [Ba/Fe] as large as 0.7. However, very high [Ba/Fe] 
(> 0.9) at [Fe/H]^ —0.6 can not be explained by the 
present starburst models. These stars with very high 
[Ba/Fe] could have been formed directly from AGB ejecta 
that did not mix well with the ISM. 

3.2.4. fg-[Fe/H] relation 

Figure 10 clearly shows that the burst model with the 
Salpeter IMF can not reproduce simultaneously the ob- 
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served fgfi and [Fe/H]o owing to cfBcient chemical enrich- 
ment. On the other hand, a larger amount of gas needs 
to be consumed for the ISM to be chemically enriched to 
[Fe/H]~ -0.3 in the models with steep IMFs (a = 2.75) 
so that /g.o can become significantly smaller in the model. 
The models with a = 2.55 yet different star formation 
histories in the quiescent phases can best reproduce both 
the observed /g,o and [Fe/H]o. These results on the IMF- 
dependences are essentially the same as those derived for 
the standard model. 

3.3. Double-burst models 
3.3.1. AMR 

Figure 11 shows the dependences of the AMRs on IMFs 
and star formation histories in the quiescent phase for 
the double-burst models. The derived IMF dependence 

of AMRs arc very similar to those in the burst ones. The 
model DB3 with steep IMF (a = 2.75) and less rapid star 
formation in the quiescent phase is inconsistent with any 
observational results on the AMR. The model DBl with 
the Salpeter IMF can better reproduce the observed AMRs 
of C05 and COS whereas the model DB2 with moderately 
steep IMF (a = 2.55) can be better fit to the AMR by H09. 
Although the burst at t = 12.8 Gyr is as strong as that 
at i = 11 Gyr, the signature of the burst in the AMR is 
less clear in the double-burst models. The observed AMR 
shows a large [Fe/H] dispersion in the youngest stellar pop- 
ulation (i.e., [Fe/H] ranging from ~ —0.6 to ~ 0), and 
the presence of such relatively metal-poor ([Fe/H]~ —0.6) 
stars at the present time is puzzling. The recent accretion 
of metal-poor gas and the resultant active star formation 
(before the burst at t = 12.8 Gyr) could result in the 
formation of such metal-poor stars and thus introduce a 
larger scatter in [Fe/H] at the young populations. 

3.3.2. [Mg/FeJ 

Figure 12 shows that the time evolution of [Mg/Fe] 
in the double-burst models is characterized by two oc- 
currences of sharp [Mg/Fe] increase/decrease after star- 
bursts, which result in two bumps. The double-burst mod- 
els show two bumps in the [Mg/Fe]-[Fe/H] plane, though 
such bumps are not clearly seen in the observations. The 
final [Mg/Fe] is larger than for all double-burst models 
owing to the last starburst at t = 12.8 Gyr, which is a 
feature that discriminates between the burst and double- 
burst models. The models (DB3, 4 and 5) with steeper 
IMFs show smaller final [Mg/Fe] and the final [Mg/Fe] 
in the model DBl with a = 2.35 is more consistent with 
the observed one for the LMC field stars. The observed 
[Mg/Fe] of the field stars with -0.5 <[Fe/H]< -0.3 ap- 
pears to increase from ~ to 0.2, which may be regarded 
as the second bump due to the starburst at t = 12.8 Gyr. 
The second bump can be shown more clearly in the models 
DBS and DBS with steeper IMFs (a = 2.75). As noted 
in §3.2, the observed large [Mg/Fe] dispersion at a given 
[Fe/H] makes it difficult for us to confirm the presence (or 
absence) of the first and second starbursts in the observed 
[Mg/Fe]-[Fe/H] relation. 

3.3.3. [Ba/FeJ 



Figure 13 shows that [Ba/Fe] rapidly decreases soon af- 
ter the first starburst a.tt = 11 Gyr and then increases un- 
til the commencement of the second starburst a.t t = 12.8 
Gyr for the double-burst models with different IMFs (in 
the time evolution of [Ba/Fe]). The increase of [Ba/Fe] af- 
ter the second starburst can not be seen, because the 0.2 
Gyr difference between the present and the last starburst is 
not long enough for AGB stars to have chemically polluted 
the LMC ISM. Owing to the first starburst, [Ba/Fe] can 
become significantly high (up to ~ 0.7) in the double-burst 
models. However, the maximum value across all models 
is still significantly lower than those (> 0.9) observed in 
some stars with [Fe/H]= —0.6 ^ —0.5. The models with 
steeper IMFs can show higher final [Ba/Fe] and they can 
better explain the clear differences in the locations of field 
stars in the [Ba/Fe]-[Fe/H] plane between the LMC and 
the Galaxy. 

3.3.4. fg-[Fe/H] relation 

Figure 14 shows that the best model to explain both 

/g,o and [Fe/H]o simultaneously is the one with a = 2.55 
among the double-burst models (i.e., DB2 and DB4). This 
result combined with those in the standard and burst mod- 
els strongly suggests that the IMF of the LMC should be 
moderately steeper {a « 2.55) in its long-term star forma- 
tion history. As demonstrated, the models with steeper 
IMF can better explain the presence of the LMC field stars 
with [Ba/Fe] > 0.5 and higher [Ba/Fe] of GCs at low metal- 
hcities ([Fe/H]< —1.5). It should be noted here that all 
cjccta from AGB stars and SNe are retained in the LMC 
for these non-wind models. Thus, it can be conchided that 
a steeper IMF (a ~ 2.55) can better explain the star for- 
mation and chemical evolution histories of the LMC, as 
long as the LMC has not lost a significant amount of its 
chemically enriched ISM through stellar winds. 

3.4. Wind models 
3.4.1. AMR 

Figure 15 shows that (i) the AMRs are not so different 

between selective wind models (Wl— W4) with different 
a, Cq, and /oj and (ii) chemical evolution proceeds signif- 
icantly more slowly in the non-selective wind model (W5) 
than in the selective ones until recently (t ~ 8 Gyr) so that 
the AMR in the non-selective wind model shows system- 
atically lower [Fe/H] for a given age. Although the non- 
selective wind model can explain the observed [Fe/H]o, the 
total gas mass ejected from the LMC (Mej^t) for 13 Gyr 
is about 3.4 times larger than the final stellar mass (i.e., 
Mej,t = 9.1 X IO^Mq for = 2.7 x lO^M©) and thus 
appears to be too large. Both already existing ISM and 
newly synthesized metals can be efficiently removed from 
the LMC in the non-sclcctive wind models so that the 
chemical evolution of the LMC can proceed much more 
slowly. As a result of this, a much larger amount of gas 
can be removed from the LMC until the stellar metallic- 
ity finally becomes [Fe/H]ps 0.3 in the non-selective wind 
model (W5). The AMRs in the selective wind models are 
very similar to the standard models: They are more con- 
sistent with the observed one by C05 and C08 than that by 
HZ09. This result implies that the observed AMR alone 
does not enable us to discuss the effects of stellar winds in 
the LMC chemical evolution. 



8 



3.4.2. [Mg/Fe] 

Figure 16 shows that the monoHthic decrease of [Mg/Fc] 
with time and the [Mg/Fc]-[Fe/H] relation in the five wind 
models are essentially similar to those derived in the stan- 
dard models with no wind. The selective wind models 
with larger /oj show slightly higher [Mg/Fc] for a given 
IMF whereas those with steeper IMFs show lower [Mg/Fc] 
for a given /oj . The five wind models can not explain the 
observed higher [Mg/Fe] (> 0.2) for [Fe/H] > -0.6 in the 
LMC field stars. There is no remarkable difference in the 
[Mg/Fc] evolution between the selective and non-selective 
wind models. 

3.4.3. [Ba/Fe] 

Figure 17 shows that the selective wind models can show 
higher [Ba/Fe] (> 0.5) for [Fc/H]> —0.6 for a range of a 
and /ej. The model W2 with a = 2.55 and /ej = 0.4 shows 
[Ba/Fe] > 0.6 without secondary starbursts, because Fe- 
rich gas is ejected from the LMC disk whereas the AGB 
ejecta containing s-process elements (e.g., Ba) can be re- 
tained and thus used for chemical evolution. These re- 
sults imply that efficient removal of SNe ejecta from the 
LMC disk is a way to significantly increase [Ba/Fe] without 
starburst(s). The models with steeper IMFs show higher 
[Ba/Fe], which is essentially the same as the results de- 
rived for the standard models. The final (thus maximum) 
[Ba/Fe] in the non-selective wind model at [Fe/H] = —0.3 
is lower than 0.4, and therefore inconsistent with the ob- 
served vahie, which implies that the non-selective wind 
model is less promising than the selective wind one as a 
reasonable model for LMC chemical evolution. 

3.4.4. fg-[Fe/H] relation 

Figure 18 shows that the selective wind model Wl with 
the Salpcter IMF appears to better fit the observed /g,o 
and [Fe/II]o than the models with steeper IMF {a = 2.55) 
for /oj = 0.4 (W2). However, the model with steeper 
IMF and /oj = 0.2 (W4) can also better explain the ob- 
served two quantities than that with the Salpeter IMF and 
/ej = 0.2 (W3): it should be noted that the model W5 
can also explain /g,o but does not show the observed high 
[Ba/Fe] at [Fc/H]^ —0.6. These results mean that if we 
carefully choose the two parameters (a and /ej), the ob- 
servation can be well reproduced. The non-selective wind 
model with the Salpeter IMF can also reproduce both /g^o 
and [Fe/H]o reasonably well. We do not intend to dis- 
cuss the wind models with starburst(s), because the ef- 
fects of starbursts on the LMC chemical evolution are es- 
sentially the same as those already described in the burst 
and double-burst models. 

3.5. Comparison between the four different types of 

models 

We here briefly summarize the advantages and disad- 
vantages of the four different types of models in repro- 
ducing recent observational results (see Table 2 for a brief 
summary). The observed AMR can be consistent with 
most models with a reasonable set of model parameters. 
The models with steeper IMFs {a > 2.55) can better ex- 
plain the observed /g_o~[Fe/H] relation and [Ba/Fe] > —0.2 
at [Fe/H]> —1.5 in a self-consistent manner. The ap- 
parent bump in the observed [Mg/Fe] [Fc/H] relation at 



[Fc/H]^ —0.6 can be better reproduced in models with 

starbursts. The higher [Ba/Fe] - 0.5 at [Fe/H] 0.3 can 

be better explained by non-starburst models (i.e., standard 
ones) with steeper IMFs and starburst ones. Very high 
[Ba/Fe] (> 0.9) at [Fe/H]^ —0.6 can not be explained by 
any of the models in the present study. The significantly 
higher (> 0.2) [Ba/Fe] observed in some old, metal-poor 
GCs ([Fe/H]< —1.5) can not be explained by any of these 
models, either. 

The selective wind models with the Salpeter IMF can 
better explain observations than the non-wind models with 
the Salpeter IMF, though the selective wind models need 
to assume an apparently large /oj = 0.4 and can not accu- 
rately reproduce the [Ba/Fe] of the stars in the metal-poor 
GCs with [Fe/H]< —1.5. This leads us to suggest that a 
steeper IMF is required for explaining the different ob- 
served chemical properties of the LMC in a self-consistent 
manner. The removal of metals from SNe can also signifi- 
cantly increase the [Ba/Fe] (up to ~ 0.8) after a starburst 
if a large /ej > 0.5 is chosen in the selective wind mod- 
els. This implies that the selective removal of SNe ejecta 
could be partly responsible for the observed rather high 
[Ba/Fe] (-- 0.8) in the LMC field stars. Thus, there are 
two possible ways to explain the high [Ba/Fe]: One is star 
formation directly from Ba-rich AGB ejecta and the other 
is significantly efficient selective removal of SN ejecta. 

Thus, as summarized in Table 2, the four different types 
of models can explain most of the observed chemical prop- 
erties of the LMC reasonably well (except the unusually 
high [Ba/Fe]), as long as the model parameters are care- 
fully chosen. The observed AMR and overall trends of the 
[Mg/Fe]-[Fe/H] (and equally the [a/Fe]-[Fe/H]) and the 
[Ba/Fe]-[Fe/H] relations can give less strong constraints 
on the IMF and the star formation history of the LMC. 
Table 3 describes how the four key chemical properties of 
the LMC constrain the IMF and the presence or absence 
of past starbursts. 

4. DISCUSSION 
4.1. A steeper IMF? 

The present study has shown that the standard, burst, 
and double-burst models with steeper IMFs {a ~ 2.55) 
can explain well the observed /g and [Fe/H]o in a self- 
consistent manner. Furthermore such non-wind models 
can not only better explain the observed higher [Ba/Fe] at 
[Fe/H]< -1.5 in the LMC (due to the steeper IMFs which 
raise the r-process/Fe ratio), but also can reproduce well 
the overall dependence of [Ba/Fe] on [Fe/H]. Although the 
selective wind models with the Salpeter IMFs {a = 2.35) 
can explain both /g and [Fe/H]o in a self-consistent man- 
ner for a larger /ej (~ 0.4), they can not explain so well 
the observed [Ba/Fe] at [Fe/H]~ —1.5. Therefore, the non- 
wind models with steeper IMFs seems to be slightly better 
models for the adopted stellar yields. 

Although it seems to be premature to observationally 
determine whether the IMF of the LMC is clearly steeper 
than the Salpeter one, a number of previous observations 
have suggested a steeper IMF in the LMC. For example, 
Mateo (1988) investigated the IMFs of stars with masses 
ranging from 0.9Mq to 10.5Mq in the six clusters of the 
LMC and found that the slopes are typically a ~ 3.5. Hill 
et al. (1994) also found a steeper IMF {a 3) for young 
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stars with masses larger than 9Mq in the LMC and also 
suggested that the IMF can bo different below and above 
9Mq. Holtzman et al. (1997) investigated the IMF for 
stars on the main sequence which are fainter than the old- 
est turnoff and discussed the possibility of a steep IMF 
with a = 2.75 for the LMC dominated by young popula- 
tions. As shown in Figures 4 and 6 of this paper, rather 
steep IMFs (e.g., a > 2.95) cannot be consistent with the 
observed chemical properties and AMR of the LMC. 

The only difference in the non-wind models with steeper 
IMFs and the selective wind ones with the Salpeter ones 
is the [Fe/H]-dependence of [Ba/Fe] at lower [Fe/H] (< 
— 1.5). Most of the observed field stars and GCs with 
[Fe/H]< -1.5 in the LMC show [Ba/Fe] > -0.2, which is 
more consistent with the non-wind models with steeper 
IMFs. However the total number of data points in the 
observed [Ba/Fc] [Fe/H] diagram is sufficiently small that 
we can not make a robust conclusion as to whether the 
non-wind models or the selective wind ones are better in 
explaining the observed [Ba/Fe] trend with [Fe/H]. Thus 
it is undoubtedly worthwhile for future spectroscopic ob- 
servations to investigate [Ba/Fe] and other s-process el- 
ements of the metal-poor stellar populations of the LMC 
with [Fe/H]< —1.5 to obtain better constraints on the IMF 
and the ejection processes of SN ejecta in the LMC. 

4.2. Evidence for prompt SNe la and jet-induced SNe II? 

The present study first investigated how a prompt SNe 
la influences the chemical evolution of the LMC and 
thereby predicted that [a/Fe] decreases monotonically 
(until a secondary starburst occurs) with no remarkable 
plateau in the [Q;/Fe]-[Fe/H] relation for low [Fe/H]. We 
should first investigate [Mg/Fe]-[Fe/H] relations, firstly 
because this element is most reliably determined by the 
observations among a-elements, and secondly because its 
yield is most reliably predicted by nucleosynthesis calcula- 
tions. The second best a-element for this investigation is 
Ca. On the other hand, Ti and Si are not so good since Ti 
is not well predicted by SN II nucleosynthesis, and the ob- 
served determination of the Si abundance involves a large 
uncertainty. 

The observed [Mg/Fe]-[Fe/H] relation does not show the 
predicted trends so clearly partly because of the large dis- 
persion in [Mg/Fe] at each [Fe/H]. Although the physical 
reasons for the large [Mg/Fe] dispersion are not so clear, 
one of the posssible explanations is described as follows. 
If different lo(;al regions have different star formation his- 
tories (e.g., due to different local gas densities and gas 
infall rates) in the LMC, as observed in recent observa- 
tions (e.g., HZ09 and R12), then they can have different 
[Mg/Fe] owing to different chemical enrichment histories. 
Therefore, the observed large scatter can be due to differ- 
ent evolutions of [Mg/Fe] in different local regions of the 
LMC. 

The observed [Ca/Fc]- [Fe/H] relation can be also used 
for discussing whether there is evidence for the prompt 
SNe la playing a significant role in the chemical evolu- 
tion of the LMC. Figure 19 shows the locations of the 
LMC stars on the [Ca/Fe]-[Fe/H] plane as well as the pre- 
dicted [Ca/Fe]-[Fe/H] relations for three different models. 
These three models are chosen because they, in combina- 
tion, show widely different [Ca/Fe] [Fe/H] relations of the 



LMC stars. It is clear that [Ca/Fe] decreases as [Fe/H] in- 
creases without showing a clear plateau for [Fe/H]< —1 in 
the LMC stars, which is consistent with the predictions of 
the present models with prompt SNe la. In addition, the 
observed locations of the LMC stars on the [Ti/Fe]-[Fe/H] 
plane clearly show a trend for a wide range of metallic- 
ity similar to [Ca/Fe] (e.g., C12). However the [Si/Fe]- 
[Fe/H] relation does not show such a clear trend owing to a 
larger [Si/Fe] dispersion for stars with [Fe/H]< —1.5 (e.g., 
Figure 4 in C12). These results on the [Ca/Fe] -[Fe/H] 
and [Ti/Fe]-[Fe/H] relations are supporting evidence that 
prompt SNe la have influenced the chemical enrichment 
history of the LMC. 

Figure 19 also shows that there are significant differ- 
ences in the distribution of field stars on the [Ca/Fe]- 
[Fe/H] plane between the LMC and the Galaxy. The LMC 
stars with [Fe/H]> —1 have systematically lower [Ca/Fe] 
in comparison with their Galactic counterparts with sim- 
ilar metallicities. The significantly lower [Ca/Fe] in the 
LMC stars can not be easily explained by the models even 
with very steep IMFs with a = 2.95, though the models 
with steeper IMFs can show lower [Ca/Fe]. On the other 
hand, the observations do not show significant differences 
in the locations of field stars on the [Mg/Fe]- [Fe/H] plane 
between the LMC and the Galaxy. Given that the mod- 
els with rather steep IMFs (a > 2.75) show significantly 
low [Mg/Fe] (< —0.2), the models with unusually steep 
IMFs (q > 2.95) are miable to explain both the observed 
distributions of the LMC field stars on the [Ca/Fe]-[Fe/H] 
and [Mg/Fe]-[Fe/H] planes. So how can we explain these 
observations? 

If Ca can be expelled more effectively from the LMC 
ISM during explosions of SNe II in comparison with other 
Qf-elements, then [Ca/Fe] becomes significantly lower than 
[Mg/Fe] as chemical evolution proceeds. We here suggest 
that nucleosynthesis of jet-induced SNe II can be respon- 
sible for the origin of the rather low [Ca/Fe] as follows. 
Shigeyama et al. (2010) have recently investigated hy- 
drodynamical processes and nucleosynthesis in jet-induced 
SNe and derived aspherical distributions of chemical yields 
of the SNe. They have found that both the [Ca/Fe] and 
the [Mg/Fe] in the ejecta of a jet-induced SN II can de- 
pend strongly on azimuthal angles 9 (where = corre- 
sponds to the direction of the jet). They have also found 
that [Mg/Ca] can be lower in lower 6 in their A3 model 
in which the total explosion energy of an aspherical su- 
pernova is 10^^ erg and the chemical abundances of O, 
Mg, Fe, Ca, Cr, Mn, and Zn for each 9 are investigated 
(see their Figure 1 for the 9 dependences). If the jet- 
induced SNc II eject with lower 9 (where a larger amount 
of Ca-rich gas exists) can bo more efficiently expelled from 
the LMC, then [Ca/Fe] in the LMC becomes more rapidly 
lower (in comparison with [Mg/Fo]) as the chemical en- 
richment proceeds. Although it is not clear at this stage 
whether or not this more efficient removal of Ca from the 
LMC is really possible, we here discuss how much more 
efficiently Ca should be removed in order to explain the 
observed [Ca/Fe]-[Fe/H] relation. 

We have investigated the selective wind model (W6) in 
which Ca is removed more efficiently by a factor of 1.75 
from the LMC in comparison with other elements of SNe 
II ejecta. Therefore, My, in the equation (5) is rewritten 
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in the model W6 as follows: 

M =1 l-75/ejMej,sn for Ca 

^ \ /ojMoj^sn for other elements ^ ' 

The parameter values of a, Cgq, fej in model W6 are 
the same as those adopted in model W2 (i.e., a = 2.55, 
Cq = 0.015, and fcj = 0.4). Figure 19 demonstrates that 
model W6 shows a steeper [Ca/Fc] decrease with [Fe/H] 

and a rather low final [Ca/Fc] ( 0.3) at [Fe/H] 0.3. 

Furthermore, it is confirmed that the [Mg/Fe]-[Fe/H] in 
the model is also consistent with observations. Although 
the results of the model arc broadly consistent with obser- 
vations, it is not clear why some of the intermediate-age 
LMC GCs with [Fe/H]~ -0.5 have higher [Ca/Fe] (~ 0) 
than the field stars with similar metallicities. The possi- 
ble difference in [Ca/Fe] between the GCs and the field 
stars in the LMC could be related to the differences in 
the formation processes between field stars and GCs and 
thus is beyond the scope of this paper. We will discuss 
the origin of this intriguing difference in our forthcoming 
papers based on chemodynamical numerical simulations of 
the LMC. 

4.3. Chemical signatures for the past bursts of star 
formation 

Recent observational studies have investigated the AMR 
for different local regions in the LMC and thereby dis- 
cussed the spatially resolved star formation and chemi- 
cal evolution histories of the LMC (e.g., HZ09 and R12). 
HZ09 have found that there are peaks of star formation at 
roughly 2 Gyr, 500 Myr, 100 Myr, and 12 Myr ago in the 
LMC. R12 also have revealed the presence of peaks in the 
star formation rates at 2.0 Gyr and 250 Myr ago. Recent 
numerical simulations on the formation of the Magellanic 
Stream have shown that the LMC and the SMC could 
have experienced strong tidal interactions at about 2 Gyr 
and 250 Myr ago, and suggested that the two interactions 
could have significantly enhanced star formation in the 
LMC (DB12). These recent observational and theoreti- 
cal studies imply that the LMC might have experienced a 
burst of star formation at least twice, though the epoch 
and the strength of each burst have not been precisely de- 
termined yet. In the following discussion, we focus on the 
possible starbursts at ~ 2 Gyr and ~ 200 Myr ago. 

As shown in previous theoretical models, past starburst 
events can be imprinted on the chemical abundances of the 
stellar and gaseous components of the LMC (e.g., Russell 
& Dopita 1992; T95). In particular, [a/Fe] as a function 
of [Fe/H] can significantly change during starburst events 
(e.g., T95) and thus can be used to give strong constraints 
on the epochs of the events. PT98 clearly showed that (i) 
the non-selective wind model with a secondary starburst 
about 2 Gyr ago can better explain the observed AMR 
and (ii) [a/Fe] can significantly and rapidly increase dur- 
ing the starburst and then slowly decrease to the solar 
value (e.g., [Mg/Fc]^ 0). The present study has predicted 
that the chemical signatures of the starburst about 2 Gyr 
ago include (i) a rapid increase of [a/Fe] followed by a 
rapid decrease, (ii) a rapid decrease of [Ba/Fe] followed by 
a rapid increase, and (iii) a time delay between the [a/Fe] 
and [Ba/Fe] peaks. 

However, as pointed out in previous sections, these three 
chemical signatures of the past starburst event around 



2 Gyr ago can not be clearly seen in the observational 
results. Figure 20 shows [Mg/Fe]-[Fe/H] and [Ba/Fe]- 
[Fe/H] relations for the five burst models (B4— B9) with 
different epochs and strengths of starburst. The results 
shown in this figure (and Figure 8) suggest that although 

the observed higher [Mg/Fe] (> 0.2) at [Fe/H] 0.6 can 

be consistent with the rapid increase of [Mg/Fe] due to 
the starburst event, the large scatter in [Mg/Fe] around 
[Fe/H]~ —0.6 does not allow us to make a robust conclu- 
sion on the origin of the higher [Mg/Fe]. Similarly, the 
observed larger [Ba/Fe] (> 0.5) for [Fe/H]> -0.4, which 
is consistent with the present starburst models, can not be 
regarded as strong evidence for the presence of a starburst 
about 2 Gyr ago owing to the small number of stars with 
known [Ba/Fe] for [Fe/H] > —0.4. More observational 
data sets are necessary to investigate how the secondary 
starburst about 2 Gyr ago might have changed the chem- 
ical evolution history of the LMC. The AMRs derived by 
different observational studies are significantly different so 
that the comparison between the observed AMRs and the 
simulated ones in the present study can not provide strong 
constraints on the nature of the past starburst events ei- 
ther. 

Concerning a possible starburst about 0.2 Gyr ago, a 
chemical signature of the starburst about 0.2 Gyr ago is 
the sudden increase of [Mg/Fe] around [Fe/H]~ —0.4, as 
shown in the present double-burst models. There is a hint 
of such a [Mg/Fe] increase in the observational data for the 
LMC field stars by FOB, though the number of the stars is 
very small (only three). If the starburst is strong enough, 
then it can be imprinted on the MDF of the LMC. Fig- 
ure 21 describes the MDFs for [Fe/H] and [Mg/Fe] in the 
five double-burst models. The MDFs here are relative fre- 
quency histograms that are binned with 0.1 dex bin width 
and normalized to the most populated bin. Clearly the 
MDF of [Fe/H] show the double peaks for the five models 
and such double peaks are not observed in the previous 
observation by COS. The observed apparent lack of a bi- 
modal MDF in COS could be due to the small number of 
young and metal-rich samples in COS. If the observational 
result is real, then it means that the starburst about 0.2 
Gyr ago is much weaker than modeled in the present study 
(so that it can not be detected in the observed MDF). 

4.4. The origin of the dip in the AMR 

The AMR derived for the LMC by HZ09 (in their Figure 
20) appears to show a sudden and significant decrease of 
[Fe/H] around 1.5 Gyr ago followed by a rapid increase of 
[Fc/H], though HZ09 did not discuss the origin of the pos- 
sible "dip" in the AMR. The AMRs for some local regions 
of the LMC derived by R12 also appear to show the dips 
whereas the AMR by COS does not clearly show the dip. If 
the observed dip around 1.5 Gyr ago (HZ09) is real, it has 
profound implications on the gas accretion history of the 
LMC. One of likely explanations for the possible dip is that 
a large amount of external metal-poor gas ([Fe/H] < —1.0, 
significantly smaller than the gaseous metallicity of the 
LMC about 2 Gyr ago) was accreted onto the LMC from 
outside the LMC disk and consequently [Fe/H] rapidly and 
significantly decreased. The observed apparently rapid de- 
crease of [Fe/H] by almost 0.2 dex at [Fe/H]~ -0.7 (HZ09) 
can give some constraints on the amount of gas accreted 
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onto the LMC for a given metallicity of the metal-poor 
gas. 

In order to discuss how the AMR of the LMC can change 
owing to the infall of metal-poor gas from outside the LMC 
disk, wc have investigated the AMRs of models with gas 
infall yet without any starburst before infall. The pur- 
pose of this investigation is to show clearly how the dip 
of the AMR can be formed during the infall of metal-poor 
gas (thus not to reproduce the observed AMR fully self- 
consistently) . Therefore, wc think that the adopted some- 
what idealized models would be enough to show clearly 
the formation of the dip in the AMR due to the infall of 
metal-poor gas. We have mainly investigated how much 
gas needs to be accreted onto the LMC by using the re- 
sults of the "infall models" in which metal-poor gas with 
[Fe/H]=— 1.6 can be accreted onto the LMC at 1.5 Gyr 
ago. Here the metallicity of [Fc/H]= —1.6 is chosen just 
for a representative case to clearly show the formation of 
the dip in the AMR (we also investigated the infall models 
with [Fe/H]= —1.0 for comparison). 

In the infall models, the following infall rate of external 
gas (Mext) is added to the right side of the equation (1): 

= , ^1^1 > (12) 

''in,c ^'in,s 

where Mext is the total mass of the external gas that can 
infall onto the LMC for tin ^ < t < tin,e, where tin,s and 
tin,o are the epochs when the gas infall starts and ends, re- 
spectively. Therefore, the gas infall rate is assumed to be 
steady and constant in the present infall models. We in- 
vestigate the models with tin,s = 11-5 Gyr and tin,e = 12.0 
Gyr. The chemical abundance patterns of the gas is as- 
sumed to be the same as that of the LMC halo. Therefore, 
the term of Z/^jMext (where Zj^i is the chemical abundance 
of each element in the infalling gas) is added to the right 
hand side of equation (2) to calculate the evolution of Zi. 
Owing to the infall of metal-poor gas, the mean metallic- 
ity of the youngest population becomes significantly lower. 
Therefore, a starburst needs to occur after the gas infall so 
that the final [Fe/H] can be —0.3, as observed. We thus as- 
sume that starbursts can occur before and after accretion 
of the metal-poor gas onto the LMC disk in the infall mod- 
els, because the models with this assumption can be con- 
sistent also with the observational results by HZ09 (i.e., a 
possible starburst about 2 Gyr ago). We have particularly 
investigated the four infall models (11—14) with a = 2.55 
and Cq = 0.004, Cshi = 0.8, t,bi s = H-O Gyr, Ui c = H-l 
Gyr, U2,s = 12.0, U2,c = 13.0 Gyr, [Fc/H]= -1.6 in ex- 
ternal infalling gas, and different Csb2 and Mext- These 
models were chosen because they can together show dif- 
ferent degrees of sudden [Fe/H] drop (or different depths 
of the dip) in the AMRs. For comparison, wc have investi- 
gated a model (15) with [Fe/H]= —1.0 in external gas for 
comparison. The model parameters for these models are 
shown in Table 4. 

Figure 22 shows the AMRs for the last ^ 5 Gyr (i.e., 
<; = 8 — 13 Gyr) in the four infall models in which Mext 
ranges from 0.1 to 0.4 in simulation units. Mext = 1 thus 
means that the total amount of external gas accretion is 
the same as the total amount of gas accreted onto the 
LMC from its own halo for the last 13 Gyr. The value 
of Csb2 in each model is chosen such that the final [Fe/H] 
can be consistent with the observed one. For compari- 



son, the AMR of the burst model B2 without gas infall 
is shown in this figure. Clearly stellar [Fe/H] can rapidly 
decrease soon after the metal-poor gas is accreted onto 
the LMC with the magnitude of the decrease depending 
on the amount of the accreted gas (Mext). The model 13 
with Mext = 0.3 can show the dip of ~ 0.2 dex, which 
means that ^ 10^ Mq needs to be accreted onto the LMC 
for explaining the observed dip. Other models with lower 
Mext (i.e., II and 12) show less remarkable dips in the 
AMRs and thus are less consistent with the observations. 
If the metallicity of the accreted gas is higher than the 
adopted one ([Fe/H]=— 1.6), then even a larger amount of 
gas needs to be accreted onto the LMC to form the remark- 
able dip in the AMR. For example, if the infalling gas has 
[Fe/H]= —1, then Mext ~ 0.4 is required for explaining the 
dip of ~ 0.2 dex observed in the AMR. We thus conclude 
that if the observed dip is real, then a massive accretion 
event of metal-poor gas with Mext of at least ^ IO^'Mq is 
required to explain the observed magnitude of the dip. 

Previous numerical simulations by BC07 and DB12 
demonstrated that the gas of the SMC can be accreted 
onto the LMC after tidal stripping of the SMC gas caused 
by the strong LMC SMC Galaxy interaction. These sim- 
ulations showed that there could be two accretion events 
around 1.5 Gyr and 0.2 Gyr ago, which means that the 
first accretion event about 1.5 Gyr ago is a promising can- 
didate which can provide a large enough amount of gas 
to form the dip in the LMC AMR. However, the total 
amount of the SMC gas transferred to the LMC about 1.5 
Gyr ago in these simulations is less than IO^Mq for the 
total SMC mass of 3 x 1O^M0. The predicted gas mass is 
much smaller than the required mass (Mext ^ IO^Mq) for 
explaining the observed dip. This means that the SMC gas 
accreted onto the LMC is unlikely to explain the observed 
dip of 0.2 dex, if the total mass of the SMC about 2 Gyr 
ago is similar to that of the present SMC (^ 3 x 10^ Mq). 
However, if the SMC was originally much more massive, 
the required amount of gas mass could be transferred from 
the SMC to the LMC. 

Then where does such a large amount of metal-poor gas 
come from? One possible scenario is that a gas-rich dwarf 
galaxy with a total gas mass of IO^Mq merged with the 
LMC about 1.5 Gyr ago and the mixing of the gas with 
the LMC ISM caused a significant decrease of [Fe/H]. Since 
the gas-mass of the dwarf should be similar to that of the 
LMC in order to explain the observed dip, the dwarf would 
have to have had a total mass similar to the LMC. This 
means that the LMC stellar disk could have been severely 
damaged by the violent dynamical process of such a major 
merger event. The thick disk and counter-rotating stellar 
components (e.g., Subramaniam & Prabhu 2005) might 
have been formed from this major merger event occurring 
in the LMC about 1.5 Gyr ago. A merged dwarf as mas- 
sive as the LMC should have a much lower star formation 
rate than the LMC so that it can have a low metallicity. 
This discussion depends on the assumption that the ob- 
served dip of 0.2 dex is real and the dip was caused by 
a massive gas accretion event. Given that the amplitude 
of the dip provides information about the total mass of a 
gas-rich dwarf merging with the LMC (or the total mass of 
gas accretion onto the LMC), it would be quite important 
for future extensive observational studies to confirm the 
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presence or the absence of the dip(s) in the AMR. 

4.5. Formation of very high [Ba/Fe] st,ars in the LMC 

Although the present models can reproduce reasonably 
well the overall trend of [Ba/Fe] with [Fe/H] and the higher 
[Ba/Fe] (- 0.5) for [Fe/H]< -0.6 in the LMC, they can 
not explain the stars with [Ba/Fe] > 0.9. We have con- 
firmed that even the selective wind models with steep 
IMFs and efficient metal ejection (/ej 0.6) can show 
at most [Ba/Fe] ~ 0.8. This failure of the present one-zone 
chemical evolution models is related to the adopted as- 
sumption that AGB cjccta can be mixed with ISM soon 
after AGB stars eject their s-process elements. If the AGB 
ejecta does not mix well with the surrounding ISM and 
consequently can be converted into new stars, then the 
stars can show rather high [Ba/Fe]: Figure 1 in TB12 
shows the observed rather high (> 1) [Ba/Fc] in the en- 
velope of AGB stars with different metallicities. We thus 
propose that the stars with unusually high [Ba/Fe] (> 0.9) 
in the LMC were formed as a result of incomplete mixing 
of ISM and AGB ejecta. A key question here is how en- 
ergetic stellar winds from AGB stars can cool down to 
become cold gas for star formation without mixing so well 
with the surrounding ISM. 

Recent hydrodynamical simulations have shown that 
secondary star formation directly from AGB ejecta is pos- 
sible in massive star clusters owing to the deeper gravi- 
tational potentials (e.g., Bekki 2011). Although this sec- 
ondary star formation appears to be a convincing mecha- 
nism for the formation of stars with unusually high [Ba/Fc] 
in the LMC GCs, it (;a,n not explain why some of the LMC 
field stars show such high [Ba/Fe]. One possibility for the 
field star formation from AGB ejecta is that AGB ejecta 
can assemble in the HI holes (e.g., Kim et al. 1999), where 
ISM can be almost completely blown away by SNe, and 
then can be converted into cold gas there without mixing 
efficiently with chemically enriched ISM. The new stars 
formed from AGB ejecta in the HI holes can naturally 
have very high [Ba/Fe]. Thus it would be important for 
observational studies to confirm whether the locations of 
young stars with very high [Ba/Fe] are more likely to be 
in the present HI holes. Our future chcmodynamical sim- 
ulations will investigate whether this star formation from 
AGB ejecta in HI holes is really possible in the LMC. 

5. CONCLUSIONS 

We have investigated the chemical evolution of the LMC 
by using new one-zone chemical evolution models in which 
both chemical pollution by prompt SNe la and metallic- 
ity dependent chemical yields of AGB stars are incorpo- 
rated and the IMF is a free parameter. We have particu- 
larly investigated three different types of models with or 
without starburst (i.e., standard, burst, and double-burst 
models) so that we can discuss the importance of previous 
starbursts in the history of the chemical evolution of the 
LMC. We furthermore have investigated the wind models 
in which gaseous ejecta from AGB stars, SNe la, and SNe 
II can be removed partly from the LMC owing to stel- 
lar feedback effects of SNe and AGB stars. The principle 
results of the models are summarized as follows. 

(1) The observed gas mass fraction (/g^o ~ 0.3) and 
the metallicity of youngest stellar populations ([Fe/H]o ~ 



—0.3) in the LMC together give some constraints on the 
IMF and the efficiency of gas removal by stellar feedback 
effects. Both /g and [Fe/H]o can be best reproduced by 
a steeper IMF with a = 2.55 for the standard (i.e., no 
burst) models in a self-consistent manner, if the gaseous 
ejecta from AGB stars and SNe are not removed from the 
LMC (i.e., non-wind models). This tendency of steeper 
IMFs to better explain /g and [Fe/H]o simultaneously can 
be seen also in the burst and double-burst models. Fur- 
thermore, the observed higher [Ba/Fe] (> —0.2) of stars in 
GCs at [Fe/H]< —1.5 in the LMC is also more consistent 
with steeper IMFs (a > 2.55). 

(2) However, the models with the Salpetcr IMF (a = 
2.35) can also explain the observed /g^o and [Fe/H]o, if 
significant fractions (/ej ~ 0.4) of gaseous ejecta from SNe 
are selectively removed from the LMC (i.e., no removal of 
AGB ejecta) . These selective wind models are significantly 
more reasonable than the uniform wind ones proposed by 
PT98 in which ISM, AGB ejecta, and supernova ones are 
equally removed. This is because an unreasonably large 
amount of gas (~ lO^^M©, or three times the present stel- 
lar mass of the LMC) need to be removed from the LMC 
for the best uniform wind models with the Salpcter IMF. 
The wind models with the Salpeter IMF, however, can 
not reproduce well the observed higher [Ba/Fe] (> —0.2) 
at lower [Fe/H] (< -1.5). Thus, although removal of SNe 
ejecta could be important in the chemical evolution of the 
LMC, an IMF steeper than the Salpetcr one is required 
for explaining /g^oi [Fe/H]o, and [Ba/Fe] at low [Fe/H] in 
a self-consistent manner. 

(3) The present models predict that [a/Fe] starts to de- 
crease monotonically only ~ 10* yr after the commence- 
ment of active star formation in the LMC owing to rapid 
chemical enrichment by ejecta from prompt SNe la. The 
observed [Ca/Fe]-[Fe/H] relation with the apparent lack 
of the plateau can be consistent with models with chem- 
ical pollution by prompt SNe la rather than by classical 
ones. This suggests that the observed [Ca/Fe]-[Fe/H] re- 
lation is supporting evidence for prompt SNe la playing a 
key role in the chemical evolution of the LMC. However, 
it should be noted that [Mg/Fc] docs not so show such a 
clear trend as decreasing [Mg/Fe] with increasing [Fe/H] 
(like [Ca/Fe]-[Fe/H] relation) for the entire sample of stars 
and GCs. 

(4) The present models predict that if the LMC ex- 
perienced a starbiirst about 2 Gyr ago, [a/Fe] can start 
to rapidly increase (up to 0.3) at [Fe/H]~ —0.5 owing to 
gaseous ejecta of SNe II and then soon decrease owing to 
the Fe-rich ejecta from prompt SNe la. Therefore, the best 
model predicts a "bump" in the [q!/Fc]-[Fc/H] relation 
around [Fe/H]~ —0.5, though such a bump can not be so 
clearly seen in the observed [a/Fe] [Fe/H] relation owing 
to the apparently large scatter of [a/Fc] at [Fe/H]~ —0.5. 
The observed presence of stars with [Mg/Fe] > 0.3 at 
—0.6 <[Fe/H]< —0.5 and the observed apparent lack of 
stars with [Mg/Fe]< 0.2 at -0.5 <[Fe/H]< -0.3 can be 
consistent with the presence of the bump, though the 
smaller number of observational data points for [Fe/H]> 
—0.5 could be responsible for the apparent dip. 

(5) If the LMC experiences a secondary starburst, 
[Ba/Fe] temporarily decreases owing to Fe-rich ejecta from 
SNe II and prompt SNe la then increases sharply owing to 



13 



AGB ejecta that is rich in s-proccss elements. For exam- 
ple, if the LMC experieneed a starburst about 2 Gyr ago 

when the LMC had [Fe/H] 0.7, then the LMC shows 

the [Ba/Fe] peak at [Fe/H] 0.3 in the [Ba/Fc]-[Fc/H] 

relation. Furthermore, [Ba/Fc] can show its peak (up to 
~ 0.7) always after [a/Fe] takes its peak value. As a natu- 
ral result of this, [Fe/H] at the [Ba/Fe] peak is higher than 
that at the [a/Fe] peak in the [Ba/Fe]-[Fe/H] and [a/Fe]- 
[Fe/H] relations. These predicted trends of [Ba/Fe] are 
not so clearly seen in the observed [Ba/Fe] -[Fe/H] rela- 
tion owing to the large scatter of [Ba/Fe]. However, the 
observed larger [Ba/Fe] at [Fe/H] ~ —0.3 is more con- 
sistent with the present burst models (in particular with 
those of steeper IMFs). However, the observed very large 
[Ba/Fe] ~ 1 at [Fe/H]~ —0.5 in some field stars and GCs 
of the LMC can not be simply explained by any model in 
the present study. Thus we have proposed that such very 
high [Ba/Fe] stars could be formed from AGB ejecta that 
did not mix well with ISM. 

(6) The observed AMR has a large scatter, so the stan- 
dard (non-burst), burst, and double-burst models can be 
all consistent with the AMR, which implies that the AMR 
does not give strong constraints on the LMC star forma- 
tion history. The present double-burst models have bi- 
modal distributions of [Fe/H] (i.e., two peaks in the MDF), 
which have not been observed yet. Therefore, the double- 
burst models are the least consistent with observations 
among the; thrcic different types of models investigatcid in 
the present study in terms of the MDF. Accordingly, we 
suggest that the observationally inferred recent starburst 
around 0.1-0.5 Gyr ago by HZ09 should be weak so as to 
reproduce the observed MDF. 

(7) If the observed apparent dip (i.e., a sudden [Fe/H] 
decrease by ~ 0.2 dex) in the AMR around 1.5 Gyr ago 
is real, then it has a profound implication for the gas 
accretion history of the LMC. The present models pre- 
dict that the dip could be due to the accretion of a large 
amount (~ IO^Mq) of metal-poor gas ([Fe/H]< —1) from 
other gas-rich galaxies about 1.5 Gyr ago. Given that 
previous numerical simulations (e.g., BC07 and DB12) 
demonstrated a gas transfer from the SMC to the LMC 
caused by tidal stripping of the SMC about 1.5 Gyr ago, 
the SMC gas could be responsible for the gas accretion 
event in the LMC. However, the required large amount of 
gas (~ lO^M©) is much larger than the predicted amount 



(■^ 10^ Mq) in previous simulations. There can be two pos- 
sible scenarios for the gas transfer. One is that the SMC 
was originally much more massive than the present SMC, 
as suggested by the recent modeling of the SMC's rotation 
curve (Bekki & Stanivirovic 2009), so that a large amount 
of gas can be transferred from the SMC to the LMC. The 
other is that the LMC experienced a major merger with 
a massive gas-rich dwarf galaxy. Such a gas-rich major 
merger event may be responsible for the formation of the 
thicker and extended stellar disk observed in the LMC. 

(8) The observed rather low [Ca/Fe] (< -0.2) at 
[Fe/H]> —0.6 in the LMC field stars can not be simply 
explained by the present models, in which all elements of 
SNe are equally efficiently removed from the LMC. We 
have thus proposed that if Ca can be (by a factor of ^ 2) 
more efficiently removed from the LMC through supernova 
feedback effects in comparison with other elements, then 
the observed unusually low [Ca/Fe] and normal [Mg/Fe] 
for [Fe/H]> —0.6 can be simultaneously explained. Al- 
though it remain unclear why such more efficient removal 
of Ca should occur in the chemical enrichment history of 
the LMC, we have proposed that a characteristic nucle- 
osynthesis of jet-induced SNe can be associated with the 
origin of the field stars with unusually low [Ca/Fe]. The 
origin of the higher [Ca/Fe] in the intermediate-age GCs 
of the LMC remains unclear. 

(9) The differences in IMFs and removal efficiencies of 
AGB and SN ejecta between the LMC and the Galaxy 
might be responsible for the observed clear differences in 
the locations of the stars in the [a/Fe]-[Fe/H] and [Ba/Fe]- 
[Fc/H] relations between the two galaxies. Our future 
chemodynamical simulations of the star formation and 
chemical enrichment histories in the LMC will investigate 
how and why the IMF and the removal processes of stellar 
ejecta of the two galaxies might be different. We plan to 
discuss the spatially different chemical properties in the 
LMC based on the results of the simulations. 
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Table 1 



Model parameters for one-zone chemical evolution 



Model 


a ^ 






^sbl,s 


^sbl,e 


Csb2 ^ 


isb2,s ^ 


tsb2,e ' 


J ej ^ej 


SI 


2.35 


0.006 


- 


- 


- 


- 


- 


- 


- 


S2 


2.15 


0.004 


- 


- 


- 


- 


- 


- 


- 


S3 


2.55 


0.01 


- 


- 


- 


- 


- 


- 


- 


S4 


2.75 


0.017 


- 


- 


- 


- 


- 


- 


- 


S5 


2.95 


0.02 


- 


- 


- 


- 


- 


- 


- 


Bl 


2.35 


0.004 


0.3 


11.0 


11.1 


- 


- 


- 


- 


B2 


2.55 


0.004 


0.8 


11.0 


11.1 


- 


- 


- 


- 


B3 


2.75 


0.004 


1.2 


11.0 


11.1 


- 


- 


- 


- 


B4 


2.55 


0.006 


0.5 


11.0 


11.1 


- 


- 


- 


- 


B5 


2.75 


0.006 


1.4 


11.0 


11.1 


- 


- 


- 


- 


B6 


2.35 


0.004 


0.03 


11.0 


12.0 










B7 


2.55 


0.004 


0.8 


10.0 


10.1 










B8 


2.55 


0.004 


0.8 


8.0 


8.1 










B9 


2.35 


0.004 


0.6 


6.0 


6.1 










DBl 


2.35 


0.004 


0.2 


11.0 


11.1 


0.2 


12.8 


12.9 




DB2 


2.55 


0.004 


0.55 


11.0 


11.1 


0.55 


12.8 


12.9 




DB3 


2.75 


0.004 


1.0 


11.0 


11.1 


1.0 


12.8 


12.9 




DB4 


2.55 


0.006 


0.4 


11.0 


11.1 


0.4 


12.8 


12.9 




DBS 


2.75 


0.006 


0.9 


11.0 


11.1 


0.9 


12.8 


12.9 




DB6 


2.35 


0.004 


0.15 


11.0 


11.1 


0.15 


12.5 


12.6 




Wl 


2.35 


0.01 














0.4 


W2 


2.55 


0.015 














0.4 


W3 


2.35 


0.008 














0.2 


W4 


2.55 


0.012 














0.2 


W5 


2.35 


0.012 














0.2 300.0 


W6 


2.55 


0.015 














0.4 



'^Thc "S", "B", "DB", and "W" are referred to as the standard (i.e., no burst), burst, double- 
burst, and wind models, respectively. The removal efficiency /ej is by a factor of 1.75 higher in 
Ca than other a-elements in the wind model W6. 

^The slope of the IMF. 

'^The coefficient for star formation in the quiescent (i.e., no starburst) phase. 

•^The coefficient for star formation in the first starburst phase. 

•^The time at which the first starburst begins in units of Gyr. 

^The time at which the first starburst ends in units of Gyr. 

^The coefficient for star formation in the second starburst phase. 

^The time at which the second starburst begins in units of Gyr. 
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' The time at which the second starburst ends in units of Gyr. 

j The mass fraction of SN ejecta that can be completely removed from the LMC in the selective wind models. 
^ The coefficient for the metal ejection rate in the non-selective wind models. 



Table 2 



Comparison of different models 



Properties/models 


Standard 


Burst 


Double-burst 


Wind 


AMR 





o 


o 





[Mg/Fe] - [Fe/H] 








o 





a bump in the [Mg/Fe] — [Fe/H] relation 


X 





o 


0^ 


[Mg/Fe] at [Fe/H]o 





o 


o 





[Ba/Fe] - [Fc/H] 





o 


o 





high (- 0.5) [Ba/Fc] at [Fc/H] - -0.3 


0^ 


0^ 






very high (> 0.9) [Ba/Fe] at [Fe/H] ~ -0.6 


X 


X 


X 


X 


[Ba/Fe](> -0.2) of GCs at [Fe/H] < -1.5 










/g,o - [Fe/H]o 








o 



*If a listed chemical property of the LMC is (not) reproduced reasonably well by a model, 
then a diagnosis mark "Q" ("x") is given for the property. 

''A starburst about 2 Gyr ago is required for the (selective and non-selective) wind models 

'^The IMF slope a needs to be steeper than 2.55 in these standard, burst, and double-burst 
models. 

■^Only selective wind models can reproduce this observation. 

*^Thc IMF slope a needs to be steeper than 2.55 in these four different types of models. 
^The IMF slope a needs to be ~ 2.55 in these standard, burst, and double-burst models. 
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Table 3 

coxsth \ints ok the imf am) the stak fokmatiox lustoin' l\ the laic 

Properties * Requirement 

a bump in the [Mg/Fe] — [Fe/H] relation starburst about 2 Gyr ago 

high 0.5) [Ba/Fe] at [Fc/H] 0.3 a > 2.55 in non-wind models ^ 

[Ba/Fe](> -0.2) of GCs at [Fe/H] < -1.5 a > 2.55 

/g,o ~ [Fs/HJq a ~ 2.55 in non-wind models 



''The listed four properties are selected from the seven in Table 2 because 
they can give some constraints on the IMF and the star formation history of 
the LMC. 

''This property can be reproduced in wind models if only SNe ejecta are 
partly removed from the LMC (i.e., only in selective wind models) for a rea- 
sonable range of IMFs. 



Table 4 



Parameters for the infall models 



Model 


Mext 


[Fe/H] in external gas 


C'sb2 


11 


0.1 


-1.6 


0.05 


12 


0.2 


-1.6 


0.07 


13 


0.3 


-1.6 


0.1 


14 


0.4 


-1.6 


0.12 


15 


0.4 


-1.0 


0.12 
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Time (xio* yr) 



Fig. 1. — The time evolution of star formation rates (SFRs) in units of Mq yr~^ (upper) and [Fe/H] (lower) for the five representative 
models v^ith different star formation histories, SI (blue solid), Bl (cyan dotted), B6 (red short-dashed), DBl (magenta long-dashed), and 
DB 6 (green dot-dashed). The SFRs are estimated for these models by assuming that the present total stellar mass of the LMC (Mg) is 
2.7 X IO^Mq. 
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Fig. 2. — Chemical evolution of four models on the [Mg/Fe]-[Fe/H] plane (upper) and the time evolution of [Mg/Fe] for two prompt SN la 
models with /i, = 0.03 (blue solid) and 0.06 (cyan dotted) and two classical SN la ones with = 0.03 (red short-dashed) and 0.06 (magenta 
long-dashed). 
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Fig. 3. — The AMRs for the five standard models (S1-S5) with different a (IMF slope): a = 2.15 (blue solid), a = 2.35 (cyan dotted), 
a = 2.55 (red short-dashed), a = 2.75 (magenta long-dashed), and a = 2.95 (green dot-dashed). Observational results of the LMC field stars 
by C05 (blue triangles), COS (black filled circles), and HZ09 (red pentagons) are shown. The AMR for young clusters and GCs from HZ09 
are also plotted by open circles. 
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Fig. 4. — Chemical evolution on the LMC disk on the [Mg/Fc]-[Fe/H] plane (upper) and the time evolution of [Mg/Fe] (lower) for the five 
standard models (SI— S5) with a = 2.15 (blue solid), a = 2.35 (cyan dotted), a = 2.55 (red short-dashed), a = 2.75 (magenta long-dashed), 
and a = 2.95 (green dot-dashed). The observed locations of the LMC field stars (big filled circles) and clusters (big open circles) and the 
Galactic field stars (small dots) on the [Mg/Fe]-[Fc/H] plane are shown for comparison. The observational results include P08 for the LMC 
field stars, Johnson et al. 2006 (J06), MucciareUi et al. (2008, 2010, 2011), and C12 for the LMC clusters, Gratton et al. (1999) Reddy et al. 
(2003) Venn et al. (2004), and Bensby et al. (2005) for the Galactic field stars. 
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Fig. 5. — The same as Figure 4 but for the [Ba/Fe]-[Fe/H] relations and the [Ba/Fe] evolution. 
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Fig. 6. — The same as Figure 4 but for the /g— [Fe/H] relations and the /g evolution. The filled circle with a vertical error bar indicates 
the observed present gas mass fraction of the LMC (/g,o). 
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Fig. 7. — The same as Figure 3 but for the five burst models, Bl with a = 2.35 and Cq = 0.004 (blue soHd), B2 with a = 2.55 and 
Cq = 0.004 (cyan dotted), B3 with a = 2.75 and Cq = 0.004 (red short-dashed), B4 with a = 2.55 and Cq = 0.006 (magenta long-dashed), 
and B5 with a = 2.75 and Cq = 0.006 (green dot-dashed). The parameter Cq can control the rapidity of star formation in the quiescent 
phase before starburst (see the main text for more details of Cq). 
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Fig. 8. — The same as Figure 4 but for the five burst models, Bl with a = 2.35 and Cq = 0.004 (blue solid), B2 with a = 2.55 and 
Cq = 0.004 (cyan dotted), B3 with a = 2.75 and Cq = 0.004(red short-dashed), B4 with a = 2.55 and Cq = 0.006 (magenta long-dashed), 
and B5 with a = 2.75 and Cq = 0.006 (green dot-dashed). 
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Fig. 9. — The same as Figure 5 but for the five burst models, Bl with a = 2.35 and Cq = 0.004 (blue solid), B2 with a = 2.55 and 
Cq = 0.004 (cyan dotted), B3 with a = 2.75 and Cq = 0.004(red short-dashed), B4 with a = 2.55 and Cq = 0.006 (magenta long-dashed), 
and B5 with a = 2.75 and Cq = 0.006 (green dot-dashed). 
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Fig. 10. — The same as Figure 6 but for the five burst models, Bl with a = 2.35 and Cq = 0.004 (blue solid), B2 with a = 2.55 and 
Cq = 0.004 (cyan dotted), B3 with a = 2.75 and Cq = 0.004(red short-dashed), B4 with a = 2.55 and Cq = 0.006 (magenta long-dashed), 
and B5 with a = 2.75 and Cq = 0.006 (green dot-dashed). 
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Fig. 11. — The same as Figure 3 but for the five double-burst models, DBl with a = 2.35 and Cq = 0.004 (blue solid), DB2 with 
a = 2.55 and Cq = 0.004 (cyan dotted), DB3 with a = 2.75 and Cq = 0.004(red short-dashed), DB4 with a = 2.55 and Cq = 0.006 (magenta 
long-dashed), and DB5 with a = 2.75 and Cq = 0.006 (green dot-dashed). The strength of the starbursts (Cgbi and Csb2) a^t 2 Gyr ago and 
0.2 Gyr ago (tgbi s ^ * ^ isbl,c and tsb2,s t < tsb2,c) at each model are chosen such that the final [Fe/H] is the same as the observed one 
(i.e., —0.3). The parameters tgbi s and tsbl,c (isb2,s and tsb2,c) can control the start and the end of the first (second) starburst, respectively, 
in the LMC. 
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Fig. 12. — The same as Figure 4 but for the five double-burst models, DBl with a = 2.35 and Cq 
a = 2.55 and Cq = 0.004 (cyan dotted), DBS with a = 2.75 and Cq = 0.004(red short-dashed), DB4 with a 
long-dashed), and DB5 with a = 2.75 and Cq = 0.006 (green dot-dashed). 
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Fig. 13. — The same as Figure 5 but for the five double-burst models, DBl with a = 2.35 and Cq = 0.004 (blue solid), DB2 with 
a = 2.55 and Cq = 0.004 (cyan dotted), DB3 with a = 2.75 and Cq = 0.004(red short-dashed), DB4 with a = 2.55 and Cq = 0.006 (magenta 
long-dashed), and DB5 with a = 2.75 and Cq = 0.006 (green dot-dashed). 
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Fig. 14. — The same as Figure 6 but for the five double-burst models, DBl with a = 2.35 and Cq 
a = 2.55 and Cq = 0.004 (cyan dotted), DBS with a = 2.75 and Cq = 0.004(red short-dashed), DB4 with a 
long-dashed), and DB5 with a = 2.75 and Cq = 0.006 (green dot-dashed). 
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Fig. 15. — The same as Figure 3 but for the five wind models, Wl with a = 2.35 and /cj = 0.4 (blue solid), W2 with a = 2.55 and 
/ej = 0.4 (cyan dotted), W3 with a = 2.35 and /pj = 0.2(red short-dashed), W4 with ol = 2.55 and /ej = 0.2 (magenta long-dashed), and W5 
with OL = 2.35 and Ccj = 300 (green dot-dashed). CI— C4 are selective wind models whereas C5 is a non-selective wind model. Starbursts do 
not occur in these wind models. 
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Fig. 16. — The same as Figure 4 but for the five wind models, Wl with a = 2.35 and /cj = 0.4 (blue solid), W2 with a = 2.55 and 
/ej = 0.4 (cyan dotted), W3 with a = 2.35 and f^j = 0.2(red short-dashed), W4 with a = 2.55 and f^j = 0.2 (magenta long-dashed), and W5 
with a = 2.35 and C^j = 300 (green dot-dashed). 
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Fig. 17. — The same as Figure 5 but for the five wind models, Wl with a = 2.35 and /cj = 0.4 (blue solid), W2 with a = 2.55 and 
/ej = 0.4 (cyan dotted), W3 with a = 2.35 and f^j = 0.2(red short-dashed), W4 with a = 2.55 and f^j = 0.2 (magenta long-dashed), and W5 
with a = 2.35 and C^j = 300 (green dot-dashed). 
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Fig. 18. — The same as Figure 6 but for the five wind models, Wl with a = 2.35 and /cj = 0.4 (blue solid), W2 with a = 2.55 and 
/ej = 0.4 (cyan dotted), W3 with a = 2.35 and f^j = 0.2(red short-dashed), W4 with a = 2.55 and f^j = 0.2 (magenta long-dashed), and W5 
with a = 2.35 and C^j = 300 (green dot-dashed). 
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Fig. 19. — Chemical evolution on the [Ca/Fe]-[Fe/H] plane for the three models S3 (blue solid), S5 (cyan dotted), and W6 (red dashed). 
The observational plots are from observational studies shown in Figure 4. The wind model W6 is different from other wind models in that 
Ca can be by a factor of 1.7 more efficiently removed from the LMC disk in comparison with other elements from SNe. This selective removal 
of Ca is based on previous theoretical studies on nucleosynthesis of jet-induced SNe (Shigeyama et al. 2010). Note that the wind model W6 
can be better fit to the observed low [Ca/Fe] (< -0.2) for the LMC field stars with [Fe/H]> -1.0. 
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Fig. 20. — Chemical evolution on the [Mg/Fe]-[Fe/H] plane (upper) and the [Ba/Fe]-[Fe/H] one (lower) for the five starburst models with 
different epoch and duration of starburst; B4 (blue solid), B5 (cyan dotted), B7 (red short-dashed), B8 (magenta long-dashed), and B9 (green 
dot-dashed). 
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Fig. 21.— The MDFs for [Fe/H] (left) and [Mg/Fe] (right) in the five double-burst models, DBl (top), DB2 (second from the top), DB3 
(middle), DB4 (second from the bottom), and DBS (bottom). The MDFs are normalized to the maximum numbers of stars in the bins. 
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Fig. 22. — The AMRs for the last ~ 5 Gyr for the burst model B2 (blue solid) and the four infall models, II with Afcxt = 0.1 (cyan dotted), 
12 with Mcxt = 0.2 (red short-dashed), 13 with Afcxt = 0.3 (magenta long-dashed), and 14 with A/oxt = 0.4 (green dot-dashed). Here Afcxt 
represents the total mass of external metal-poor gas accreted onto the LMC. Mext = 1 thus means that the total amount of external gas 
accretion is the same as the total amount of gas accreted onto the LMC from its own halo for the last 13 Gyr. Infall of external metal-poor 
gas onto the LMC disk is assumed to commence 1.5 Gyr ago and end 1 Gyr ago in these four infall models. Starbursts are assumed to occur 
twice at 2 Gyr ago and 1 Gyr (just after the end of the external gas infall event) in these infall models so that the final stellar metallicity in 
each model can be the same as the observed one ([Fe/H]r^ —0.3). 



